Non-existence of the final first integral in the Zipoy-Voorhees space-time 
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We show that the geodesic motion in the Zipoy-Voorhees space-time is not Liouville integrable in 
that there does not exist an additional first integral, meromorphic in the phase space variables. 



I. INTRODUCTION 

The question of integrability of the test particle mo- 
tion in the Zipoy-Voorhees metric has recently attracted 
some attention, with both numerical T, 2} and analytical 
investigations [3| . The significance of this problem is the 
same as in the classical result by Carter Q that there 
exists an additional first integral in the Kerr space-time, 
making the problem completely integrable. Carter's in- 
tegral is not generated by a Killing vector, so it is not 
a usual symmetry of the manifold, but it is quadratic 
in momenta which has important consequences. Such 
integrals translate into the separability of the Hamilton- 
Jacobi equation and d' Alembertian Q , which in turn ap- 
pears in the Teukolsky equation. That is to say, both 
the classical problem of particle motion in this space-time, 
the linear perturbation equations governing the gravita- 
tional waves and potentially quantum equations in that 
background become considerably easier to solve. This 
fact is also used in numerical approaches, when trying to 
determine possible spectra of gravitational radiation in 
anticipation of the observed data Q ■ 

It is then natural to analyse other space-times which 
could serve as models of compact objects, and the sta- 
tionary axisymmetric ones are one direction to explore. 
However, despite some numerical evidence |l[ we find 
that the particular Zipoy-Voorhees metric with the pa- 
rameter 6 = 2 is not integrable. To be more precise, we 
consider the motion of a test particle as a Hamiltonian 
system, and ask for existence of an additional constant 
of motion that would yield Liouvillian integrability with 
respect to the canonical Poisson bracket {}. That is, for 
all first integrals Ik we would have {Ik, Ii) — 0, where the 
Hamiltonian is included as H — I\. It turns out, that no 
such first integral can be found in the class of meromor- 
phic functions, and we have to use the differential Galois 
theory to prove that. 



* |mac iejka@astr o. ia.uz.zgora.pl| 

t M . Przybylska@proton. if . uz . zgora. pi 
t stachowiak@cft.edu.pl 



II. FORMULATION OF THE PROBLEM 

The Zipoy-Voorhees metric under consideration is 
aven by 



ds 2 



x - 1 



dt 2 



(x + 1) 3 (1 ~y 2 ) AA2 
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x + l 

(x 2 - l) 2 (x + l) 4 



(x 2 — y 2 ) 3 
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dx 2 
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x 2 — 1 1 — y 
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where x, y and <j) form the prolate spheroidal coordinates. 
Instead of working directly with the geodesic equations 
we take the Hamiltonian approach with 



2 y 2(>-l) 2 ^ 
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(x 2 -y 2 ) 3 (l-y 2 )^ 
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where the canonical coordinates are 



q = t, qi = x, q 2 = y, q 3 
The equations then read 

dqi dH 

dr dpi ' 
dpi dH 
dr dqi ' 



(2) 



(3) 



(4) 



and the normalization of four-velocities gives the value 
of the (conserved) Hamiltonian to be H = —\n 2 - The 
new time parameter is the rescaled proper time fir = s, 
which allows us to include the zero geodesies for photons 
without introducing another affine parameter but with 
simply /i = 0. 

Since the metric has two Killing vector fields d t and 
c?^, the two momenta po and ps are conserved. Together 
with the Hamiltonian they provide three first integrals. 
The question then is, whether there exists one more first 
integral that would make the system Liouville integrable. 
The precise notion is expressed in the following theorem. 
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Theorem 1 If the Hamiltonian system (j4]) is Liouville 
integrable then the identity component of the differential 
Galois group of the normal variational equations along 
any non-constant particular solution of is Abelian. 

Let us try to explain the involved mathematics some- 
what, although for a proper, detailed exposition of the 
subject the reader is referred to [8l4Io| , and a short intro- 
duction with application to another relativistic system 
can be found in [11| ■ 

The key steps are that we can find a particular solution 
ip (r) of our equations, and then construct the variational 
equations around it. Writing (|4]) jointly as u\ = fi(u), 
we define the variational equations as 



i= — 

duj 



(5) 



The normal part is the subspace formed by excluding 
the variations tangent to the particular solution and pro- 
jected onto the constant energy hypersurface, which re- 
duces the number of equations by 2. In general, one often 
changes the independent variable r — > z for practical rea- 
sons so that the normal variational equations (NVE) can 
be written as 



$(z)=A i i (z)£ J (z). 



(6) 



Since the first integrals are to be meromorphic as func- 
tions of positions and momenta, note that we are now 
dealing with complex functions, and extend also the time 
to the complex domain. So even though the particular so- 
lutions and equations might lie outside the physical (real) 
domain, the non-existence of complex first integrals will 
preclude the existence of real ones as well, since the latter 
are just special cases of the former. 

The integrability of the original system translates into 
solvability of the variational equations, which are much 
easier to analyse being linear. It is at this level, that 
we define the group of transformations acting on some 
function field to which the solutions of the NVE belong. 
The subgroup which leaves the of the elements of A? (z) 
fixed is the differential Galois group Q. 

Theorem 2 // the differential Galois group of a linear 
system of equations has a solvable element of identity, 
then all the solutions of the system are Liouvillian. 

In other words, if the group in question has an Abelian, 
hence solvable, identity component, the NVE are solv- 
able with Liouvillian functions. Fortunately, one does 
not have to check the group directly, and in the case of 
second order equations with rational coefficients there is 
an effective algorithm for finding the Liouvillian solutions 
.12]. 

The plan of attack is thus to look for particular solu- 
tions, for which the NVE has a block structure so that a 
two dimensional subsystem can be separated. We then 
rewrite is a second order linear differential equation with 



rational coefficients and apply the Kovacic algorithm to 
see if it has any Liouvillian solutions. Note that the 
system has no external parameters, and only the values 
of particular first integrals enter as internal parameters. 
They are synonymous with initial conditions, so that if 
we manage to find just one solution, for particular values 
of /i, po and P3, such that the respective NVE is unsolv- 
able, we will have proven that there cannot exist another 
first integral over the whole phase space. 

It might so happen, unlike in the Carter case, that the 
system exhibits some particular invariant set on which 
there exists an additional integral. For example one could 
have J4 = H, which would mean that 1 4 is conserved on 
the zero energy hypersurface fj, — which is clearly a 
physically distinguished case. We will then have to look 
for particular solutions on those sets to make the results 
even more restrictive than just the lack of a global first 
integral. 



III. PROOF OF NON-INTEGRABILITY 

The obvious particular solution to look at is a particle 
moving along a straight line, through the center in the 
equatorial plane, which in prolate coordinates means y = 
and P3 — 0. The non-trivial equations then read: 



ch _ (x + l) 2 p Q 
dr " 
d.r 



(x- l) 2 
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Pi a; 



dr (x-l)0 + l)' c 
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dpi 2 a; 5 (3-2a;) 2 2(x + l) 



Pi- 



dr rx (x + l) 6 {x - l) 2 Po \x- l) 3 ' 
Or, upon rescaling time by 

dr = r du, 



(8) 



we have 



PiX 3 (x — 1) 
(x + 1) 2 ' 



Pi =Pv 



2 x 2 (3-2x) 2 2(a; + l) 4 



(9) 



{x + lf Pa 'x 3 (x-1) 



where the dot denotes differentiation with respect to u, 
and we have omitted the t equation, as the other two 
separate from it. This is the two-dimensional subsystem 
of the NVE that was mentioned before. The conservation 
of the Hamiltonian now reads 

1 2 _ p 2 (x + lf+p 2 x & (l-x 2 ) 

2 M ~ 2(x-l) 2 (x + l) e ' 1 j 

which together with the equation for x yields 

x 2 = (x 2 - l)(p 2 (x + l) 2 - M 2 (x - l) 2 ), (11) 



3 



so that x(u) is expressible by the Jacobi elliptic functions. 
This fact is important, as we will change the independent 
variable from u to x which is permissible (does not change 
the identity component of the Galois group) only if the 
function x{u) defines a finite cover of the complex plane 

!. 

The variational equations along this solution separate 
so that the NVE read 



In case 2, the families E c are 



6 



(1+4 



• _ 2 x(x-l) 

(x + \y 

and can be brought to the standard form of 

R(x) 



(12) 



" 4x 2 (x 2 - l) 2 (p 2 (x + l) 2 - n*{x - l) 2 ) 2 
=: r{x)£ >2 {x), 



Ux) 



(13) 



where R is the following polynomial 

R{x) =Po(34x 2 - 40a; + 3) (a: + l) 4 

- 6p^ 2 (6x 2 - 10a; + l)(.x 2 - l) 2 (14) 
+ ^{22x 2 - 2Qx + i){x - l) 4 . 

Since for all physical particles we have po ^ all the 
others parameters can be rescaled by it (fx — > fx/po, p3 — > 
pB/po), which we use in what follows. 

A. General r(x) 



The poles of r(x) are 



-1,0,1 



p, — 1 fi + l 

II + 1 ' // — 1 J ' 



(15) 



and for all of them to be different we must have /i 2 ^ 1 
and n 7^ 0. All are of order 2, and the order at infinity is 
4, so that we need to check all the cases of the Kovacic 
algorithm. 

In case 1, the appropriate quantities form the fol- 
lowing set 

f(0 l) (l -§.) (l (3 1\ (5 _1\ (5 

\V u :- L yiV4' 4/'\2' 2/'V4'4/'V4' 4/'V4' 4//' 

(16) 

(the first pair corresponding to oo), and the combinations 



oo / . ^ c 

C.3 



(17) 



give only nine non-negative integers (not all distinct) as 
possible degrees of the appropriate polynomial P. How- 
ever the respective test solutions require that [i = 0, and 
have to be discarded so that this case cannot hold. 



{(0,2, 4), (9, 2, -5), (6, 2, -2), 
(3, 2,1), (5, 2,-1), (5, 2,-1)}, 



(18) 



which in turn give 131 possible integer degrees for the 
polynomial P. Checking them one by one, we find that 
they require fj, = ±1 in order to form a solution, so that 
this case can be discarded as well under the current as- 
sumptions. 

In case 3, the families E c contain 6 x 13 = 78 numbers, 
which make 4826809 combinations for d out of which 
230856 are non-negative integers. We thus first, resort to 
checking for the presence of logarithms in the solutions, 
which would prevent this case [9j . 

The only poles with integer difference in the exponents 
are and oo. Using the Frobenius method, we get the 
two independent solutions around zero 



Vi =x 



3/2 h 4- 5 ^ 2 ~ 2 ) + 23// 4 - 38// 2 + 65 2 
3(1 -M 2 ) 12(1 -/i 2 ) 2 X 

l? - 1 5(^ 2 - 2) 



v-2 = ln(a;)«i + x 



9(5 -At 2 ) 9(5 -a« 2 )' 



(19) 



As can be seen, they require that p? ^ 5, and when the 
equality holds the coefficients in the second series have 
to be recalculated to give 



Vl = x 



v 2 = x 



■ 1 1 1,. . 7 A, 



-1/2 



4 32 

, 15 1485 
1 + — x- 
4 



64 



(20) 



The solutions around oo do not have logarithms at all, 
which means that the only possibility for case 3 left here 
is with fi 2 = 5. 



B. Special subcases 

In order to exclude the special energy hypersurfaces 
we have to resort to a more general particular solution, 
namely one with p% ^ 0. As already mentioned, it is 
enough to find one solution for each such surface, and 
that means we can take a specific value of p%. The cor- 
responding NVE will only have numeric coefficients, and 
checking for its Liouvillian solutions is much easier, for it 
suffices to use the "kovacicsols" routine implemented in 
Maple []. 

The solution will also be expressible by (hyper-)elliptic 
function as defined by the Hamiltonian constraint 



{x-l){x + lf -{x-lfp 2 -{x 2 -If 



(x + l) 2 



(21) 
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and the counterpart of the NVE (|12|) will read 



£1 = 

e 2 = 



(1 + .x) 3 
(x - 1) {Zp\x A 



(x 2 - l) 2 p§) 



(22) 



£ 3 (:r + l) 2 



a, 



We then proceed exactly as above, taking x as the inde- 
pendent variable and reducing the system to one equation 
of the form £ 2 ' = r £- F° r each hypersurface in question, 
the value of p$ = 1 leads to a NVE that is not solvable 
with Liouvillian functions. This finishes the proof for all 
possible levels of the Hamiltonian. 

To further illustrate the complexity of this system, 
we have also obtained a Poincare section for the plane 
y = shown in Fig. [TJ The numerical integrator was 
based on the Bulirsch-Stoer modified midpoint scheme 
with Richardson extrapolation. We note that the special 
solution (TTT|) lies entirely in the plane y = and is a 
trajectory beginning and ending at the singularity so it 
does not contribute to the section. It also lies outside the 
visible chaotic region, which is confined to a very small 
subset of the phase space, as mentioned in 0- 



so the right hand side is not necessarily a polynomial or 
rational function. Accordingly, the rationalization of the 
NVE might not preserve the identity component of the 
differential Galois group. However, when S is rational we 
can still proceed by taking a new dependent variable to 
be 



1 



x - 1' 



(25) 



as this leads to the normal form (|13j) which involves only 
integral powers of w and w s . Assuming then that 5 = p/q, 
we can make the NVE rational by taking w x l q as the 
new variable if need be. Unfortunately, the number of 
poles (and their values) now depends on p and q, so the 
Kovacic algorithm has to be applied to each S separately, 
but in each of them it is as straightforward as above to 
use the Maple package, once suitable numeric values of 
the parameters have been chosen. For example, we have 
verified that 6 = 1/2 is also non-integrable, confirming 
the numerical evidence of [2( that both 5 > 1 and S < 1 
do not admit additional first integrals. 



IV. GENERAL METRIC 

The above results carry, to some extent, to the general 
Zipoy-Voorhees metric given by 

ds2 = -(^) 5 d^ + (|±i)^(^- 1)( i-,W 

(23) 

The main problem that arises for arbitrary S is that the 
special solution might no longer be a (hyper-)elliptic func- 
tion, because the Hamiltonian now gives 

i 2 =1 (1 - ^) " (* + ir 25 ((* + 1) 2 V - i)pl 

-(.-^-(^-ij'+V), 

(ML 



V. CONCLUSIONS 

Our main result can be stated as 



Theorem 3 There does not exist an additional, mero- 
morphic first integral of the geodesic motion in the Zipoy- 
Voorhees metric ([T]). 

This confirms the previous considerations of 3j, and 
goes much further than first integrals polynomial in mo- 
menta. Meromorphic functions include not only the an- 
alytic functions of both momenta and coordinates, but 
also rational and transcendental ones as long as their sin- 
gularities are just poles. This result thus strongly reduces 
the possibility of using constants of motion expansion in 
solving the equations of geodesic motion or gravitational 
waves. 
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